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Abstract 

The Bethe-Salpeter formahsm in the instantaneous approximation for the interaction 
kernel entering into the Bethe-Salpeter equation represents a reasonable framework for 
the description of bound states within relativistic quantum field theory. In contrast to 
its further simplifications (like, for instance, the so-called reduced Salpeter equation), 
it allows also the consideration of bound states composed of "light" constituents. Every 
eigenvalue equation with solutions in some linear space may be (approximately) solved 
by conversion into an equivalent matrix eigenvalue problem. We demonstrate that the 
matrices arising in these representations of the instantaneous Bethe-Salpeter equation 
may be found, at least for a wide class of interactions, in an entirely algebraic manner. 
The advantages of having the involved matrices explicitly, i.e., not "contaminated" by 
errors induced by numerical computations, at one's disposal are obvious: problems like, 
for instance, questions of the stability of eigenvalues may be analyzed more rigorously; 
furthermore, for small matrix sizes the eigenvalues may even be calculated analytically. 
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1 Introduction 

The Bethe-Salpeter formalism JI], is generally accepted to represent the appropriate 
framework for the description of bound states within relativistic quantum field theory. 
Within this formalism, a bound state is described by its "Bethe-Salpeter amplitude," 
which is defined as the (time-ordered) product of the field operators of the bound-state 
constituents between the vacuum and the bound state. In principle, this bound-state 
amplitude should be found as solution of the (homogeneous) Bethe-Salpeter equation. 

However, apart from very few special cases — like the famous Wick-Cutkosky model 
which describes the interaction of two scalar particles by exchange of a massless scalar 
particle — the Bethe-Salpeter equation turns out to be practically not tractable. One of 
the main reasons for this fact is the appearance of timelike variables in the equation of 
motion. Consequently people usually consider some three-dimensional reduction of the 
Bethe-Salpeter equation. The most popular among these three-dimensional reductions 
is based on the assumption that the interaction between the bound-state constituents 
is instantaneous in the center-of-momentum frame of the bound state; the result of this 
is called the "Salpeter equation" or "instantaneous Bethe-Salpeter equation" [^ . This 
equation may be formulated as eigenvalue problem for the mass M of the bound state. 
Its dynamical quantity is the "Salpeter amplitude," obtained from the Bethe-Salpeter 
amplitude by equating the time variables of the involved bound-state constituents. 

By expanding the Salpeter amplitude into some convenient set of basis matrices in 
Dirac space [0], Lagae 0, ^ managed to reduce the full instantaneous Bethe-Salpeter 
equation for fermion-antifermion bound states to a set of coupled equations for radial 
wave functions. His easy-to-handle formalism immediately stimulated renewed interest 
in the full Salpeter equation as a manageable tool for describing mesonic bound states 
05 0; H) B Uni • The first question to be answered in this context refers to the Lorentz 
structure of the confining Bethe-Salpeter kernel. The most important outcome of these 
studies is the following: 

• Lorentz-scalar confinement turns out (1) to face serious problems with respect to 



the stability of the resulting solutions 0, 0, [TT[ and (2) to be in conflict with the 



phenomenologically dictated linearity of the mesonic Regge trajectories P, [7|, pi]] . 

• Time-component Lorentz- vector confinement, on the other hand, does not share 
either of these two problems, that is, it yields (1) stable solutions H, ||, |TT] and (2) 
linear Regge trajectories for the bound states of quark-ant iquark systems [0, pJ]] . 

The square of the norm of the amplitudes obtained as solutions of the instantaneous 
Bethe-Salpeter equation is not positive definite; in order to stay in the physical sector 
of the theory, in all these analyses one has to retain only bound states of positive norm 
squared but to reject all bound states of zero or negative norm squared. In this sector, 
all mass eigenvalues M of the instantaneous Bethe-Salpeter equation are guaranteed 
to be real |, §. 

A standard technique for determining the solutions of a given eigenvalue equation is 
the conversion of this equation into the equivalent matrix eigenvalue problem. (For the 
instantaneous Bethe-Salpeter equation, this route has been followed in Refs. 0, |^, |Tl|.) 
We show here that, for a suitable choice of basis states in the space of solutions [0, the 
matrix representation of the instantaneous Bethe-Salpeter equation may be explicitly 
given in analytical form. Among others, this might prove to be advantageous for future 
discussions of the stability of solutions of the instantaneous Bethe-Salpeter equation. 



2 The Instantaneous Bethe— Salpeter Equation 

It goes without saying that we would hke to demonstrate our procedure for solving the 
instantaneous Bethe-Salpeter equation at the simplest conceivable example. First of 
all, let us agree on two simplifying assumptions: 

• The (full) fermion propagators Si{pi, mt), i = 1,2, entering in the Bethe-Salpeter 
equation may be approximated by the corresponding free propagators, given by 

S^fiiPi, m^) = -i ilt^Pi -mi) , 1 = 1,2 . 

In the above free-propagator approximation, the masses mi, i = 1,2, of the two 
fermionic bound-state constituents are then interpreted as some effective ones. 

• The two bound-state constituents have equal masses m, that is, mi = 1712 = m. 
(The generalization to the case of unequal masses of the bound-state constituents 
is, of course, straightforward |^.) 

Moreover, we illustrate the developed technique by investigating mesons carrying the 



pion's quantum numbers, treated as bound states of a quark-antiquark pair [jl3[ . Thus, 
we set up the instantaneous Bethe-Salpeter equation for ^Sq bound states of massless 
constituents with confining interaction kernels of time-component Lorentz- vector type. 

2.1 Pseudoscalar Bound States of Massless Constituents 

We consider fermion-antifermion bound states of total spin J, with parity and charge- 
conjugation quantum numbers P and C, expressed in terms of J, given by P = (— l)"'"'"^ 
and C = (— l)"^, respectively. In usual spectroscopic notation, these states are denoted 
by ^Jj. The general expansion of the Salpeter amplitude x in terms of a complete set of 
Dirac matrices performed in Ref. |^ involves precisely eight independent components. 
For the bound states under consideration, only two of these independent components 
are relevant. We call the corresponding wave functions in momentum space \E'i and "^2- 
In the notation of Lagae ||^, these wave functions are labelled Li and L2, respectively. 
For two particles of equal masses m and internal momentum k constituting a bound 
state with the quantum numbers specified above, the Salpeter amplitude x describing 
this bound state in its center-of-momentum frame reads in momentum space 



X(k) 
with the abbreviation 



^i(k)^^^^^^ + v&2(k)7° 
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E{k) = Vk^ + m2 , A;= |k| . 
The norm ||x|| of the Salpeter amplitude x may be calculated from Eq. (2.9) of Ref. 



[or from Eq. (9) of Ref. [11 1). For the above amplitude, the squared norm is given by 
in accordance with Eq. (4.13) of Ref. M, and with Eq. (18) of Ref. |TTI 



'^"' = V tS^ ^'^'^^'^ '^'^^'^ + '^^^'"^ '^'^^'^^ 



For simplicity, we focus our interest to quark-antiquark bound states of spin J = 0; 
these are bound states with the spin-parity-charge conjugation assignment J^^ = 0~~^. 
It is these pseudoscalar quark-antiquark bound states which will be investigated here. 



In view of the findings recalled in the Introduction, let us settle on the investigation 
of the instantaneous Bethe-Salpeter equation for some time-component Lorentz vector 
interaction between the bound-state constituents, that is, on a Bethe-Salpeter kernel 
of the Dirac structure j^^j^. For this case, the radial wave functions \l/i(A;) and '$2{k), 
obtained from \l/i(k) and \l/2(k) by factorizing off the spherical harmonics involving the 
angular variables, satisfy the following system of equations — which is equivalent to the 
instantaneous Bethe-Salpeter equation and which may be deduced, e.g., from Eq. (5.7) 
or Eq. (5.9) of Ref. §, or from Eq. (Al) of Ref. ^■. 



2E{k)^2{k) + 



' U'2 



dk'k 



Voik,k')^2ik') = M^,{k) 



2E{k)'^i{k) 



(1) 



+ 



/ U'2 



dk'k 



Vo{k,k) 



E{k) 



h h' 

+ ^V,ik,k') 



E{k') E{k) 



E{k') 



^i{k') =M^2{k) 



where, in terms of some interaction potential V{r) in configuration space, 

00 
VL{k,k') = 8TT fdrr^V{r)jL{kr)jL{k'r) , L = 0, 1, 2, . . . . (2) 



Here, jn{z) {n = 0, ±1, ±2, . . .) are the spherical Bessel functions of the first kind |I^ . 
For a vanishing mass of the two bound-state constituents (that is, for m = 0), this 
set of equations simplifies to 

00 

2k^il2ik) + —-^ldk'k'^Voik,k')^2ik') = M^H.ik), 





00 



2k^i{k) + j^Jdk'k'^V,{k,k')^i{k') = M^2{k). 



(3) 



The above set of equations constitutes the starting point for the present investigation. 



2.2 Analytical Solution by Series Expansion 

The particular structure of the instantaneous Bethe-Salpeter equation (|l]) allows us to 
extract its solutions by the following simple procedure. Express one of the components 
\l/i, \l/2 of the Salpeter amplitude x from one of the coupled equations in the system ([^) 
by the other component. Insert this expression into the other equation of this system, 
in order to obtain an eigenvalue equation for the other component with eigenvalue M^: 

M'^^2ik) = 4A;2^2(fc) 



+ 2k 



I Ui2 



dk'k 
72^ 



V,{k,k')^2{k') 



+ 2 



' ^'3 



dk'k 



Vi{k,k')^2{k') 



+ 



' Wi 



dk'k 



Vi{k,k') 



II vn 



dk" k 



K.(fc',fc")^2(n- 



Clearly, all our results should be independent of which component has been selected in 
the first place. For the present analysis, we have chosen to express \E'i in terms of \E'2- 
An approximate solution to the above eigenvalue equation may be very easily found 
by expanding the "Salpeter components" over a suitable set of basis functions H, |6|, pi]] . 
To this end, we introduce complete ortho normal systems {\(j)i),i = 0,1,2, ...} as bases 
for the Hilbert space L2{R^) of (with the weight function w{x) = x^) square- integrable 
functions f{x) on the positive real line R^; these basis vectors \(f)i) satisfy {(f)i\(j)j) = 6ij. 
The different sets of basis vectors are distinguished by some parameter £ = 0,1,2,.... 
The representations of any basis vector |0j) in configuration and momentum space are 
related to the corresponding bases for the Hilbert space L^^B?) of all square- integrable 
functions on the three-dimensional space E? by factorizing off the spherical harmonics 
y^miS^) for angular momentum £ = 0,1,2,... and its projection m = —i, —i+l, . . . ,+i 
which describe the dependence on the angular variables summarized by the solid angle 
Q. The parameter i is identified with the angular momentum of the R^ basis functions. 
For a given value of i^ the R^ basis functions in configuration and momentum space are 
called (pi (r) and 0^ (p), respectively. They satisfy the orthonormalization conditions 

oo 

drr'<pf\r)<pf\r) = Jdpp'<pf\p)<pf\p)=S,,, z, j = 0, 1, 2, . . . . 



In order to introduce an additional degree of freedom in the search for the solutions 
of the instantaneous Bethe-Salpeter equation, we allow the basis functions to depend 
on some positive real variational parameter n (with the dimension of mass). Our choice 
of radial basis functions, involving, in their configuration-space representation (pl (r), 
the generalized Laguerre polynomials [0], is summarized in due detail in Appendix ^ 
A particular feature of this choice is our phase convention: all configuration-space basis 
functions (pl (r) and, for i even, the momentum-space basis functions (p] (p) are real. 
Starting with expansions over the radial basis functions (pl (p) corresponding to i = 0, 
the solution of the instantaneous Bethe-Salpeter equation then simply amounts to the 
diagonalization of the — a priori, infinite-dimensional — matrix (with eigenvalues M^) 

Mij = Aij + Bij + dj + Di^ , (4) 

with the abbreviations 

oo 

A,^ ^ 4 f dk k' <pf\k) <pf\k) , (5) 



Bij 



°° °° 1 , / 7/2 

2Jdkk'ci\k)j^V,{k,k')<pf{k') , (6) 



(27r) 

^ ^ 



oo oo 



dk' k'^ 



a, ^ 2jdkk'<pf\k)J—--V,{k,k')<pf\k'), (7) 



(27r) 

^ ^ 



~'V;./2 '^AU"U'2 



A, = Jdkk'<pr\k)J^V,{k,k')J^^Vo{k',k")<pf\k"). (8) 

In actual calculations, an infinite-dimensional matrix, arising from a countably infinite 
number of basis states \4>i), is not manageable. One has to be content with a truncation 
of the matrix Aiij to a, say, dx d matrix. Moreover, in the expansions encountered in 
the intermediate steps only the first, say, A^+1 basis vectors can be taken into account. 



Tlie (straightforward) evaluation of the matrix Aiij is briefly sketched in Appendix |B|: 

N N N 

A^,, = 4 ij^ (^) + 2 ^ 6„ (/x) kT (/") + 2 E E <. ds, (/i) Ki'^ (/i) 



r=0 

N N N 

+ EEE<.c.,v^«(/i)v^,f(/x). 

r=0 s=0 i=0 



r=0 s=0 



Let US present the various quantities entering in this result in their order of appearance. 
The (real and symmetric) matrix j|j- (/i) may be easily read off from the result (^) for 

the integral l\j{^) defined by Eq. (^) for the case n = 2: 

4;u2 



l\f{^^) 



TT ^{i + 1) (z + 2) (j + 1) (j + 2) 



EE(-2) 



r+s 



i + 2 



r=0 s=0 

■|r-s| / 

fc=0 \ 



r+s+4 

- E 

fc=0 



J + 2 



|r — s\ 
k 



(r + l)(s + l) 
r(i {k + 3)) r(i (r + s + 1 + |r - s| - A;)) 



kn 



r(|(r + s + 4+ |r-s|)) 



■ cos 



r + s + 4 \ r(|(A; + 3))r(|(2r + 2g + 5-A;)) ^^^ /^tt^ 
A; ) r(r + s + 4) ^° V 2 ^ 



Note the scaling behaviour of this expression with respect to the variational parameter 

j|)(/i) = /iV,?Hl). 



Explicitly, the matrix /'•^-'(/i) = f/j^ (/i)j reads 

/ 



I^{^^)^{Ij^\^^))=^^ 



2 
73 



2 

73 

7 

3 



V 



/ 



According to Eq. (^), the real and symmetric matrix bij{^) is identical to the integral 

ij 



Hj{p) defined by Eq. ( ^2]) for the case n = 1; consequently, it may be read off from the 



suit (53) for 


n = l: 


hij{ii) = - 


^ix 


Tv'(z + l)(^ + 2)(j + l)(j + 2) 


X >;>;(-2)^^^ +' ^.+' (r+i)(.+i) 




\\'^\ ( r s\ \ Ta{k + 2))Ta{r + s + 2+\r-s ~ k)) fkn 




_to\ ^ J r(i(r + . + 4+ r-s|)) ^"" V 2 


1 


•^4 / + s + 4 \ r(|(A; + 2)) r(i (2 r + 2 s + 6 - A;)) /A; 7r\ 1 
t^o\ ^ 1 Tir + s + A) '°^l2JJ 



Note the scaling behaviour of this expression with respect to the variational parameter 

bij{fi) = fibij{l) . 
Explicitly, the matrix 6(/i) = {hij{n)) reads 



6(/i) = {hji^i)) 



8/i 
3^ 



1 7 
Ts 5 



\ 



J 



According to Eq. ( P3| ) , the (purely imaginary) matrix Qj is nothing else but the integral 



j\j{lj) defined by Eq. (0) specified to the case n = 0; consequently, it may be read off 
from the result (^) for n = (here, because of the singularity of r(0), some care has to 



be taken when performing the limit n ^ in the algebraic expression ( P^D for J^ (/i)): 

8 



TT ^(z + 1) (z + 2) {t + 3) (2 + 4) (j + 1) (j + 2) 



r=0 s=0 



i + 4 
i — r 



J + 2 



r + 2 



+ E 



fc=2 



-^(i (4 + r + s + |r - s|)) + ^(4 + r + s) 

r-s\\ T{lk)T{^ {4 + r + s + \r - s\ - k)) 
k J r(i(4 + r + s + |r-s|)) 



kn 



cos 



4+!:+^ / 4 + r + s\ r(|fc)r(|(8 + 2r + 2g-A:)) /fcvr ^ 
^ I A; J r(4 + r + s) ""^ A 2 , 



fc=2 
1 

r + 3 



-^(| (5 + r + s+|l + r- s|)) + ^(5 + r + s) 



^'1^'W \i+r-s\\ Tak)Ta{5 + r + s+\l + r-s\-k)) (kn' 

2^ , I „,^ ,. ^TT COS 



fc=2 



k 



r(i(5 + r + s+|l + r-s|)) 



^+1+^ ( 5 + r + s \ r(|A;)r(|(10 + 2r + 2s-A;)) (kn' 
- 2^ I , I — trr ^ cos 

fc=2 \ 



k 



r(5 + r + s) 



where '^{z) denotes the logarithmic derivative of the gamma function, the "digamma 
function." Note that the quantities q^ are independent of the variational parameter /x. 
Explicitly, the matrix c = (qj) reads 



c = C; 



16 
3^ 



Vs 
1 



1 

15 
19 



\ 



5 VIE 15^5 



Similarly, according to Eq. (Q), the (purely imaginary) matrix dij^jj) is identical to the 



integral J^J' {^) defined by Eq. (^4]) specified to the case n = 1; consequently, it may be 



read off from the result 
dijifi) 



for n = 1: 



8/i 



TT ^{i + !)(« + 2) (t + 3) (z + 4) (j + 1) (j + 2) 
X E E (-2)'^' (r + 1) (r + 2) (r + 3) (s + 1) h + ^ 



r=0 s=0 
1 



I — r 
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X 



r + 2 



|r'— s 

E 

fc=0 



r — s 



k 
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r(|(4 + r + s + |r-s|)) 



fcvr 



cos 



4^^ /4 + r + s \ r(i(A; + l))r(i(7 + 2r + 2s-A;)) (kix' 

2^ , I TTTT^ ^ COS 

fc=0 \ 



A; 



r(4 + r + s) 



1 



r + 3 
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E 

fc=0 
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fc=0 V 
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Note the scaling behaviour of this expression with respect to the variational parameter 
/i: 

dij{ii) = fidij{l) . 

Explicitly, the matrix d(fi) = (dij{fi)) reads 



d{fi) = {dij{fi)) = i- 



r{e) 



( 73 1 



1^ 



\ 



V • 



•• / 



Finally, the real and symmetric matrix V^j (/x) of expectation values of the interaction 
potential V{r) with respect to the basis states |0j) (characterized by a particular value 
of t) is defined by 



K«)/ 



V;^ (^) = {<P.\V{r)\<P,) = J drr^F(r)0r(r)0}^^(r) 





^i^)(^\A^)i 



(9) 



For power-law potentials, that is, for interaction potentials V{r) of the power-law form 



v{r) = y: 



anT 



bn 



with (arbitrary) real constants a„ and 6„, these matrix elements {(f)i\V{r)\(f)j) are easily 
worked out algebraically |T5|, |l^, ^ . The algebraic expression for the general case can 
be found in Sec. 4 of Ref. [0, Sec. 3.10 of Ref. [g, or Sec. 2.8.1 of Ref. [|17|. 

The simplest model for a confining interaction between all bound-state constituents 
is provided by a linear potential 

V{r) = \r, A>0, (10) 

whence 

oo 

VL{k, k') = 871 X J drr' J L{kr) J L{k' r) , L = 0,l,2,.... (11) 



In this case, the above-mentioned general expression for the expectation values V^j (/x) 
of the interaction potential V{r), taken with respect to our basis states |0i), simphfies 
to 



.{£) 



CI/", A) 



t\j\ 



X 



\r+s 



X 



r(2£+z + 3)r(2£ + j + 3)2/i^5 J^o 



r! si 



J 



T{2i + r + s + A) . 



W/ 



Note the scaling behaviour of the potential matrix V^^- (/z, A), £ = 0, 1, . . . , with respect 
to the variational parameter /x and the coupling strength A of the linear potential ([I0|): 



rW 



vl\^^,X) 



X 



Kf(i.i). 



The explicit expression of this potential matrix V^^\jj,, A) reads, for instance, for £ = 0, 



/ 3 



I/(°)(/.,A)^(yf(/.,A)) = A 



2/i 



-v^ 



\ 



-V3 5 



and, for 



/ 5 



V^'\^,X)^{V^\^,X)) 



2fi 



-V5 



-V^ 7 



V 



/ 



Taking into account, for the special case of a linear potential, the above-mentioned 
scaling behaviour of all the expressions entering in our main result for the matrix Aiij, 
we finally arrive at 

TV N N 



tjj 



r=0 



r=0 s=0 



\2NNN 
P r=0 s=0 t=0 



(12) 



Clearly, this structure of the matrix J^ij, with the squared masses of the bound states 
as its eigenvalues, has to emerge already on dimensional grounds: the mass dimension 
of the variational parameter // is 1 while the mass dimension of the slope A of the linear 
potential ([T0|) is 2; the potential V{r) enters into the expression (|) for the matrix Aiij 
in, at most, second order. 



3 Results 

The dependence of the matrix Aiij on the variational parameter fi, as apparent from 
Eq. (|g), is 



Mij = jj^ Eij + Fij + — Gij , (13) 

/i 



with 



E,, ^ 4/,f(l), 

TV Af TV 

F,, = 2A5:6„(l)Kf(l,l) + 2A5:$:c:,4,(l)Ki'Hl,l), 

r=0 r=0 s=0 

iV Af AT 

G., = A2EEE<.c.,\/«(i,i)v;f(i,i). 

r=0 s=0 i=0 

Any solution method relying on truncated expansions over some complete set(s) of 
basis vectors only makes sense if the approximate result obtained in this way converges 
to the exact solution of the problem under study for increasing numbers of basis vectors 
taken into account in the expansions. Assuming this convergence to hold, a reasonable 
estimate of the uncertainties induced by the truncation of the expansion series to finite 
numbers of basis vectors can be given by the inspection of the lowest-dimensional case. 
For i = j = 0, we find, for the first elements of the matrices /*^^^(1), &(1), c, and d{l), 

4?(1) = 1 , ^oo(l) = ^ > Coo = i ^ nr , c?oo(l) = i ^r ' 
and, for the expectation values V^^^{\^ 1), £ = 0, 1, of the linear potential. 

This implies, for A^ = 0, for the quantities E = Eqq, F = Fqq, and G = Gqo, 

Minimizing 

TWoo = i?E + F+—G 

/i2 

with respect to the variational parameter /i, yields 

A^oo = 2 \^EG + F . 
The above values of E, F, and G then yield, for the square of the bound-state mass M, 

M^ = ^(2 + V5)\. (14) 

For a linear interaction potential like the one used in this investigation, it would be 
natural to present the numerical results for the masses M of the bound states in terms 
of the dimensionless quantities MjyX^ as has been done, e.g., in Ref. p. However, in 
order to facilitate the comparison of our findings with the ones quoted in Ref. [11|, we 
present our results for the particular, more or less reasonable |T^ value of the slope A 
of the linear interaction potential ([T0|) adopted in Ref. [|lT|, viz., for A = 0.2 GeV^. For 
this slope we obtain, for the "one-dimensional" result (0) for the bound-state mass M , 
M = 1.69604 GeV. 



10 



In general, the approximate values of the bound-state masses M, calculated within 
the present approach as the square roots of the eigenvalues of the matrix Aiij, depend 
on the variational parameter /i. In our almost entirely algebraic method of solving the 
instantaneous Bethe-Salpeter equation, this dependence on /i has been made manifest 
in form of Eq. ([T^). However, any solution technique relying on some series expansions 
can be regarded as reasonable only if it exhibits stability with respect to the increase in 
the number of basis states, in the sense that increasing the number of basis states — in 
our case, increasing the size d of the matrix Ai^j or the number A^ of basis states taken 
into account in the series expansions performed in the intermediate steps — reduces the 
dependence of the results, the bound-state masses M, on the variational parameter fi. 

Figure |l] shows (in perfect agreement with Fig. 1 of Ref . |]lll ) the dependence of the 
masses M of the three lowest-lying bound states, i.e., the first three radial excitations, 
extracted from Eq. ([13D for a linear interaction potential (0) with slope A = 0.2 GeV , 
on the variational parameter fi. When increasing the number of basis states (by varying 
the matrix size d from d = 5 to d = 15) the formation of regions where the eigenvalues 
M become independent of // should be observed; this feature is indeed found in Fig. |l]. 
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2.5 

2 
1.5 

1 
0.5 




\L 



Figure 1: Dependence of the masses M (in units of GeV) of the three lowest-lying (positive- 
norm) J = '' bound states of two massless constituents obtained by the instantaneous 
Bethe-Salpeter equation, with a time-component Lorentz-vector confining interaction kernel 
involving a linear potential V{r) = Xr with slope A = 0.2 GeV , on the variational parameter 
/i (in units of GeV) for A^ = 49 and matrix sizes d = 5 (dashed lines) and d = 15 (full lines). 



For the ground state, the region of /i-independence corresponding to the matrix size 
d = 5 is centered somewhere around /i = 1 GeV. For this value of fi, the ground-state 
mass M, calculated for a matrix size d = 15 and A^ = 49 in the intermediate-step series 
expansions, is M = 1.656 GeV, which reveals that the rough approximation (|T^) of the 
above one-dimensional analysis is not too bad (at least for the situation studied here). 

In order to pave the way for a future quantitative discussion of the accuracy of our 
method as well as for a numerical comparison with the findings of different approaches, 
we list in Table |l| the masses M of the three lowest-lying bound states for several values 
of the variational parameter /i in the neighborhood of the above "region of stability." 
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Table 1: Eigenvalues M, in units of GeV, of the instantaneous Bethe-Salpeter equation 
for two massless (i.e., m = 0) spin-i fermions experiencing an interaction described by 
a linear potential with slope A = 0.2 GeV^ and forming bound states of radial quantum 
number n,- = 0,1,2 and spin-parity-charge conjugation assignment J^*" = 0~^ (called 
l^So, 2^So, and 3^So in the usual spectroscopic notation) as functions of the variational 
parameter /x, for matrix size d = 15 and A^ = 49 in the intermediate series expansions. 



/i [GeV 


l^So 


2^80 


S^So 


0.05 


1.807 


3.253 


5.392 


0.1 


1.673 


2.360 


3.173 


0.2 


1.656 


2.283 


2.812 


0.5 


1.656 


2.280 


2.770 


1 


1.656 


2.281 


2.771 


1.2 


1.657 


2.281 


2.771 


2 


1.658 


2.284 


2.774 


3 


1.658 


2.287 


2.843 


4 


1.662 


2.393 


3.165 


5 


1.702 


2.617 


3.614 



In our procedure, the series expansion of the radial Salpeter component '^2{k) with 
respect to the (primary) system {(pi {k),i = 0,1,2, . . .} of radial basis functions reads 



d-l 



(0), 



v^2(fc) = Ee.c^(fc). 



(15) 



j=0 



The expansion coefficients Cj in this series are obtained, in the course of diagonalization 
of the matrix Aiij introduced in Eq. (^, as the eigenvectors of Aitj. From the ffist of 
the two relations constituting the instantaneous Bethe-Salpeter equation for massless 
constituents, Eq. (§), the radial Salpeter component \l/i(A;) is then obtained in the form 



d~l 



^i{k) = J2(^^Uk) 



(16) 



j=0 



where — in order to match the structure of the ffist of Eqs. (^ — we took the liberty to 
define 



Mk) 



1 

M 



N 



0,1,2, 



2k<pr\k)+Y.n:'\i^)<pf\k) 

j=0 

The corresponding representations of the above Salpeter components "^lik) and '^2{k) 
in configuration space, denoted by \l/i(r) and "^2(1^), are calculated by a Fourier-Bessel 
transformation (for J = 0): 



^i(r) 
^2(r) 



d-l 
i=0 



d-l 



E^^C^W 



(17) 
(18) 



i=0 
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with 



V^j(r) 



1 

M 



^(0)^ 



^(0)^ 



N 
j=0 



0,1,2, 



Clearly, the relevant objects are the matrix elements calculated from these amplitudes. 
Figure ^ shows (similarly to Fig. |l| in reasonable agreement with Fig. 4 of Ref. [|l^, 
as judged by the eye) the behaviour of the radial Salpeter components \E'i(r) and '^2{f') 
in configuration space for the ground state of the physical system under consideration. 
According to Sec. p.l| , the normalization of the radial Salpeter components \E'i(A;) and 
'^2{k) in momentum space is determined by the norm ||x|| of the Salpeter amplitude x'- 



ixir 



i2ny 



dke [^\{k) ^2{k) + %{k) ^i(fc)] 



In the present approach (summarized, as far as the Salpeter amplitude x is concerned, 
in Eqs. (|16|) and (0) or Eqs. (|l^ and (0), respectively) the normalization of the radial 
component \1'2 is identical in configuration and momentum space: 



fdrr^ |^2(r)|' = fdkk^ |^2(A;)|' = Y. \^i\^ ■ 



For the plot presented in Fig. ^, we have chosen this normalization to be equal to unity. 
The normalization of the radial component \E'i then follows from Eq. (0) or Eq. (|T6|). 
Once the desired value of the norm ||x|| of the Salpeter amplitude x has been fixed, it is 
an easy task to determine the common normalization factor of \E'i and \1'2 such that the 
Salpeter amplitude involving these wave functions satisfies its normalization condition. 




Figure 2: Radial "Salpeter components" ^i(r) and ^2(^) in configuration space (in units of 
GeV^'^) for the lowest-lying (positive-norm) J^'-^ = ^ bound state of two massless spin-^ 
constituents obtained by the instantaneous Bethe-Salpeter equation, with a time-component 
Lorentz- vector confining interaction kernel involving a linear potential V{r) = Xr with slope 
A = 0.2 GeV , as a function of the relative distance r (in units of GeV~ ) of the bound-state 
constituents, for a matrix size d = 50 and A^ = 49 in the intermediate-step series expansions; 
the lower and upper curves represent the Salpeter components ^i(r) and ^2(^)1 respectively. 
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4 Summary, Conclusions, and Outlook 

The present investigation has been devoted to the exphcit construction of some matrix 
representation of the Bethe-Salpeter equation in the instantaneous approximation for 
the involved interaction kerneL As usual, this conversion to a matrix equation has been 
achieved by systematic expansions of the solutions and certain intermediate quantities 
in terms of suitably chosen sets of basis states in the Hilbert space under consideration. 
A very important feature of our particular choice of basis states is that these states may 
be represented equally well in configuration and in momentum space. As a consequence 
of this, matrix elements of operators may be evaluated in the (with respect to the given 
operator) more convenient representation space: the matrix elements of the relativistic 
kinetic-energy operator (which, because of the troublesome square root, is nonlocal in 
configuration space) can be evaluated in momentum space (where it is represented by a 
simple multiplication operator which, of course, will be reminiscent of the square root); 
on the other hand, the matrix elements of the interaction-potential operator (which is, 
in general, nonlocal in momentum space) can be evaluated in configuration space (that 
is, the environment where the interaction potentials are more frequently formulated). 
For the wide class of power-law potentials, the algebraic expression of the latter matrix 
elements can be easily deduced from the general result given in any of Refs. |jl5|, |16|, |l^. 
Merely for the sake of illustration of our technique, we adopted the linear potential as 
some kind of toy model for quark confinement. 

Our main result is the formulation of the instantaneous Bethe-Salpeter equation as 
a matrix eigenvalue problem with explicitly known matrices. We notice several highly 
welcome features of this approach, which render it a very efficient technique of solution: 



1. Because of the scaling behaviour of the involved quantities pointed out in Sec. P^ 



at every instance, all dependence on both the variational parameter /z and on the 
coupling constants entering into the interaction potential V{r) (like, in the case 
of a linear potential, the slope A) may be factored out from the matrix elements. 
As a consequence, only matrices evaluated for unit values of the parameters are 
needed in our explicit formulation of the eigenvalue problem. For a given number 
of basis states taken into account, these matrices have to be calculated only once. 

(2) 

This set of matrices includes the "moment" integrals I^j (1), all the (non-square) 
matrices of expansion coefficients &jj(l), c^, and dij{l), as well as the expectation 
values V^j (1, 1), £ = 0, 1, of the linear interaction potential (in our special case). 

2. The only numerical operation required in the course of our method of solution is 
the diagonalization of the matrix representing the instantaneous Bethe-Salpeter 
equation. For the size of this matrix representation small enough, more precisely, 
for d < 4, even this diagonalization may be performed analytically. Because of an 
adequate choice of basis, already the case d = 1 yields a remarkably good result. 

Thus, our analytic approach offers a viable way of solving the Bethe-Salpeter equation. 
Finally, we feel obliged to address the questions of the significance of the developed 
technique and to discuss prospects for an actual application of the formalism presented 
here to the description of hadrons — in particular, of mesons — as bound states of quarks 
within the quantum field theory of the strong interactions, quantum chromo dynamics. 
A quantitative discussion of these problems is certainly beyond the scope of the present 
approach. Consequently, we'll content ourselves with the addition of a few modest and 
merely qualitative comments. 
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In order to justify at all the consideration of the instantaneous approximation to 
the Bethe-Salpeter formalism, one should be able to answer the question: Which 
physical systems may be reliably described by the instantaneous Bethe-Salpeter 
equation? In other words, for which physical systems is the neglect of retardation 
effects meaningful? Any sound quantitative statement about this problem would 
require a comparison of the predictions for the bound-state masses of, on the one 
hand, the full Bethe-Salpeter equation and, on the other hand, the instantaneous 
Bethe-Salpeter equation. (This analysis would be similar to the one performed in 
Ref. [|^ for the region of validity of the reduced Salpeter equation.) However, in 
order to carry out a study of this kind, one must have at one's disposal a reliable 
method for extracting information from the full Bethe-Salpeter equation; this is, 
at present, unfortunately not the case. 

In the hopefully pedagogical example chosen for the introduction of our method, 
we confined ourselves to the (somewhat simpler) case of bound-state constituents 
with vanishing masses. Clearly, the next step in this game has to be the extension 
of this method to the case of bound-state constituents with nonvanishing masses. 
The generalization to the massive case is more or less straightforward (and can be 
done for an arbitrary Lorentz structure of the Bethe-Salpeter interaction kernel). 
Of course, the nonvanishing values of the bound-state constituents' masses entail 
a more complicated structure of the Bethe-Salpeter equation and, consequently, 
a higher degree of complexity of all resulting expressions. From our point of view, 
this fact renders the results for the general case not suitable for the introductory 
presentation of the main ideas. Accordingly, the results generalizing this study to 
massive bound-state constituents have been reserved for a separate publication, 
see Ref. |]19| ; our conclusions there are that the formalism developed here may be 



applied with due modifications to the same class of potentials as considered here. 

One long-standing problem in meson phenomenology is related to the question of 
the Lorentz structure of the confining Bethe-Salpeter interaction kernel. For the 
two reasons explained in the Introduction — namely, stability of the solutions and 
linearity of the mesonic Regge trajectories — we decided to focus our attention to 
a confining interaction kernel of time-component Lorentz-vector type. However, 
as is now well known [|, 0, |ll[], a time-component Lorentz-vector confining kernel 
leads to a (from the phenomenological point of view) wrong sign of the spin-orbit 
coupling and can thus not be considered a realistic model for the confinement of 
the colour degrees of freedom present in quantum chromodynamics. Needless to 
say, our formalism can accommodate an arbitrary Dirac structure of the kernel. 
For the most popular choices of the Lorentz nature of the Bethe-Salpeter kernel, 
in particular, for an interaction kernel of Lorentz-scalar (1 1), time- component 
Lorentz-vector (7*^(8)7*^), and full Lorentz-vector (7'^®7/j) type, the explicit form 
of the instantaneous Bethe-Salpeter equation may be read off from Refs. 0, pJ]] . 

Any realistic quark model for light mesons requires the inclusion of a mechanism 
for the spontaneous breaking of chiral symmetry. Spontaneous chiral symmetry 
breaking will manifest itself in the structure of the quark propagators; the latter 
will deviate in their form from the form relevant for the free case. In principle, in 
our formalism the inclusion of spontaneous chiral symmetry breaking in this way 
should always be possible by application of the appropriate expansions in terms 
of the basis states to the expressions involving the (non-free) quark form factors. 
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• A proper embedding of spontaneous chiral symmetry breaking should manage to 
reconcile completely the Bethe-Salpeter formalism with the Nambu-Goldstone 
theorem, which demands the existence of a massless Nambu-Goldstone boson for 
every spontaneously broken global symmetry of the theory. This would allow for 
the interpretation of all light pseudoscalar mesons vr, K, rj as (pseudo-) Goldstone 
bosons in the Bethe-Salpeter framework. (For a time-component Lorentz- vector 
Bethe-Salpeter kernel, this mechanism has been shown to operate in the desired 
way in Ref. 0. This proves that the existence of massive bound states of massless 
constituents is definitely no inherent feature of the Bethe-Salpeter formalism.) 

Thus, we may be optimistic that our technique for solving the Bethe-Salpeter equation 
in the instantaneous approximation by explicit construction of a matrix representation 
of this eigenvalue equation will eventually become an important tool for the description 
of hadrons as bound states of quarks from first principles, quantum chromodynamics. 
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A The "Generalized Laguerre" Basis 

Our choice of the basis states |0j) is defined by their configuration-space representation 

r^exp(-/ir)Lf'+')(2/xr), z = 0,l,2,..., (19) 



ki^)( 



(2/i) 



2 £+3, 



\ r(2£ + z + 3) 



(7), 



which involves the generahzed Laguerre polynomials L^ix) (for the parameter 7): the 
latter quantities are orthogonal polynomials which are defined by the power series [|14 

E(-i)' 



Lt'Hx) 



t=o 



« + 7 ] ^ 
i-t t\ 



z = 0,l,2, 



and which are orthonormalized, with the weight function x^ exp(— x), according to |n 

■'I] 1 ^5 J ^ U, i, z, . . . . 



dx x^ exp(-x) lS^^ (x) 4^^ (x) = ^^^ ^.! ^ ^^ 6, 



i\ 



The requirement of normalizability of the Hilbert-space basis states |0j) imposes on the 
variational parameter /x the constraint 

/i > . 



In this case, the configuration-space basis functions 
the orthonormalization condition 

D 

dr r^ (pi (r) (jr- (r) = 5., 



Ki^)l 



defined by Eq. (|T^), satisfy 



I] ; 



^,J =0,1,2, 



Note: the configuration-space representation of our basis states |0j) is chosen to be real. 
The momentum-space representation of these basis states, cpi (p), obtained by Fourier 



transformation of Eq. (|T9D, reads 

0f(p) ^ 



(2/i) 



2 £+3 , 



-ij p 



X 






t=0 



t! 



i + 2^ + 2 \ r(2£ + t + 3)(2^) 
i-t 



(p2_^^2)(2£+t+3)/2 



X F 



'2£ + t + 3 1+t 



p^ 



^ = 0,1,2, 



(20) 



^ 2 ' 2 ' 2'p2 + ^2 
with the hypergeometric series F, defined, in terms of the gamma function F, by '^A 

Viw 



T[u)T[v) 



E 

n=0 



T(u + n) T(v + n) z^' 



Yiw + n) 



n\ 



The momentum-space basis functions <\)\ (p) satisfy the orthonormalization condition 

'dpp'<pf\p)4\p) = S,,, z,j = 0,l,2,... . (21) 



The availability of the Fourier transform of our basis functions (pi (r) in analytic form 
represents the main advantage of our choice (plQ]). Note that the momentum-space basis 
functions are real for i = 0, as well as for all even values of t. 

^f\p)=4>f\p) for £ = 0,2,4,... . 
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B Evaluation of the Matrix M 



B.l Preliminaries 



«J 



In the course of evaluating the various matrices Aij, . . . , Dij summing up to the matrix 
Aiij according to Eq. @, we occasionally encounter ("moment") integrals of the type 



I^fifi)^ fdke^-4'\k)<pf\k) , n = 0,1,2,... . 



(22) 



An), 



According to the above definition, all the matrices I>j (/i) are both real and symmetric 
under permutation of the indices i and j. Taking advantage of a somewhat simplified 
expression for our momentum-space basis functions (p\ (p) valid in the case i = 
viz., 



M 



ip) 



4 



/i7rr(z + 3) p ^^0 

2 X -(t+2)/2 



-E(-2)*(t+l) 



X 



-^ 



sin (t + 2) arctan — 



i + 2 
i-t 

p\ 



it is rather straightforward to find the (exact) analytic expression for the integrals (|2^): 

(23) 



r(«) 



4r(/^) 



4/i" 



7r^(? + l)(2 + 2)(j + l)(j + 2) 



xEE(-2; 

r=0 s=0 



r+s 



i + 2 
i — r 



J + 2 



(r + l)(s + ll 



X 



|r-s| 

E 

fc=0 



r+s+4 

- E 

k=0 



\r-s\ \ r(| (A; + n + 1)) r(|(r + g + 3 + |r - s| - n - k)) fkji^ 
k ) r(|(r + s + 4+|r-s|)) ""^ VX 



^ + g + 4 \ r(^(fc + n + l))r(|(2r + 23 + 7-n-A;)) /^^ 
fc J r(r + s + 4) ''°'l2 



Furthermore, we will need the (exact) analytic expressions of integrals of the kind 

oo 

4f{^,)^Jdke+-<j>f\k)<pf\k), n = 0,l,2,.... (24) 



In order to get rid of the (difiicult-to-handle) hypergeometric series F entering in the 
momentum-space basis functions (j)] (p), we employ a suitable recursion formula [O 



kW 



<PY\P) 



fi^ 



X 



X 



E 

t=0 



^ vr (i + 1) (2 + 2) (i + 3) {t + 4) p2 
(-2)*/z + 4\ (t + 3)!/i* 



t! 



i-t J (p2_^^2){t+3)/2 



Vp^+1^ . ( , . ^ p\ 

sm (t + 2) arctan — 

t + 2 V /^/ ^ + 3 



sin (t + 3) arctan — 



^(1) 



Note: since (j)\ (p) is purely imaginary, all these matrices J-^ (/x) are purely imaginary. 



(n)i 



Again it is a tedious but straightforward task to work out the integrals (|2^) explicitly: 



J!f{f^ 



(25) 



'/^" 



71 ^{i + l){i + 2) (z + 3) (z + 4) (j + 1) (j + 2) 



xEE(-2rMr + l)(r + 2)(r + 3)(. + i: 

r=0 s=0 



i + 4 
i — r 



3 + 2 
3 - s 



1 



r + 2 
/ kn 



E 

fc=0 



|r - s| \ r(^ (n + A;)) r(^ (4 + r + g + |r - g| - n - A:)) 
k J r(i(4 + r + s + |r-s|)) 



X cos 



4+!:+- / 4 + ^ + s \ r(i(n + fc))r(i(8 + 2r + 2s-n-A;)) /fcvr^ 
2^ , I — ^. ^ cos 

fe=0 V 



k 



r(4 + r + s) 



1 



l+r— s 

E 

fc=0 



|1 + r — s| 
A; 



r(i (n + A;)) r(|(5 + r + s + |1 + r - s| - n - A;)) 



A;7r' 



5+r+s 

- E 

fe=0 



r(|(5 + r + s + |l + r-s|)) 
1 



cos 



5 + r + s 
A; 



rf i (n + A;)) r(i (10 + 2 r + 2 s - n - A;)) 



' kn" 



cos 



r(5 + r + s) 

The radial basis functions (f)\ (r) and (pi (p) in configuration and momentum space 
given in Eqs. ([T9|) and (pO[) , respectively, are derived from the respective basis functions 
on the (full) three-dimensional Euclidean space R^ by factorizing off the dependence on 
the angular variables. This dependence is described by the spherical harmonics yem{^) 
for angular momentum ^ = 0,1,2,... and its projection m = —i, — £+ 1, . . . , +i, which 
depend on the solid angle Q and are orthonormalized according to 

Since the two sets of basis functions on R^ are related by a Fourier transformation, the 
corresponding radial basis functions (J)] (r) and (p^ (p), i = 0,1,2, ... , are related by 



-J dpp'^ji{pr)(f)f\p) 



(26) 



di) 



ip) 



(-^y \l-^ J drr^ji{pr)(f)f\r) 





(27) 



which may be easily seen with the help of the well-known expansion of the plane waves 
in terms of spherical harmonics yim in configuration (fix) and momentum (fip) space, 

oo +£ 

exp(ip ■ x) = 47r E E ^'' 3i{pr) yLi^p) 3^fm(fix) • 

£=0 m=-i 
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B.2 Evaluation of Aij 

Expressed in terms of the integral IIj{ij) introduced in Eq. (p2D, the term Aij in the 
matrix A^jj, defined in Eq. (^, is simply given by 

B.3 Evaluation of B^j 

The term Bij in the matrix M.ij, Eq. (^), becomes, for some interaction potential V{r), 
which enters in this term linearly in form of Eq. (^ , 

oo oo oo 

Bij = - fdrr^V{r) f dk k^ (j)f\k) jo{k r) f dk' k'^ jo{k' r) (f)f {k') . 



In order to be able to apply the Fourier-Bessel relation (|26|) , we expand the expression 
kcpi \k) in terms of the momentum-space basis functions cpl (k): 

N 

k<pf\k) = Y.b,.<l>f\k). (28) 

j=0 

As a consequence of the orthonormality (^if) of the (momentum-space) basis functions 
(j)l {p), the expansion coefficients hij may be expressed in terms of the integral I^j {^) 
defined in Eq. (^: 

oo 

h, =Jdkk' # (k) <pf (k) ^ l\f (/.) . (29) 



Inserting the expansion ( pS]) involving the coefficients hij, applying the Fourier-Bessel 
relation (PB|), and remembering the definition (^ of the expectation values Vij of the 
interaction potential V{r) under consideration, the term Bij reduces to the — at least, 
for all power-law potentials algebraic — expression 

N 

Bi,=2Y.b„{f^)Kf{f^) . (30) 

T- = 

B.4 Evaluation of Cij 

The term Cij in the matrix Aiij, Eq. (^, becomes, for some interaction potential V{r), 
which enters in this term linearly in form of Eq. (^, 

oo oo oo 

dj = - fdrr^V{r) f dkk^ 4^\k) ji{kr) f dk' k'^ ji{k' r) (j)f\k') . 



In order to be able to apply the Fourier-Bessel relation (pB]), we are forced to expand 
both the momentum-space basis functions (pi (k) and the expression k (pi (k) in terms 
of the momentum-space basis functions (p] (k): 

N 

<f>f\k) = T.^n<P?ik), (31) 

j=0 

N 

k<Pf\k) = Y.d^^^?ik). (32) 
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As a consequence of the orthonormality ( pi]) of the (momentum-space) basis functions 
(pl (p), the expansion coefficients Qj and dij may be expressed in terms of the integrals 
4")(/i), n = 0, 1, defined in Eq. (11): 



c. 



«j 



.2 j,*(i)/ 7.^^(0) n.^ - 7(0) 



dA:A;Vr(A;)0f(A:) = J,;^^(/i), (33) 





00 



dij 



Jdkk'<pf\k)<pf\k)^4\f^). (34) 



Inserting the expansions (|3lD and (^) involving the coefficients Qj and (ijj, applying 
the Fourier-Bessel relation (|26|), and remembering the definition (^ of the expectation 
values V^j of the interaction potential V{r) under consideration, the term Cij reduces 
to the — at least, for all power-law potentials algebraic — expression 

N N 

Q, = 2J2T.<^ds,{^^)y|Pi^^)■ (35) 

B.5 Evaluation of Dij 

The term Dij in the matrix Aiij, Eq. (H), becomes, for some interaction potential V{r), 
which enters in this term quadratically in form of Eq. (^, 

00 00 00 

Dij = — fdrr^V{r) f dr' r'^ V{r') f dk P 0f ^ (k) ji {k r) 
00 

00 00 

X f dk' k'^ ji{k' r) jo{k' r') f dk" k"^ jo{k" r') (f)f\k") . 



^ ' I l/^\ I'm r r\-vvY\ o r\T A\^ ' i 



Adopting the expansion (RTI) of (pi (k) in terms of (pi (k) — which introduces again the 



coefficients Cij — it is clearly no problem to apply twice the Fourier-Bessel relation (|26|): 

Dij = i- ^ 4 / drr^V{r)(t)^^\r) / dr' r'^V{r') (t)f\r') / dkk'^ ji{kr) jo{kr') . 
^ ^=0 

However, in order to be able to apply the Fourier-Bessel relation (P7[), we are forced to 
expand the expressions V{r) cpl ^(r), £ = 0, 1, over the appropriate configuration-space 
basis functions cpi (r), £ = 0, 1, respectively. The corresponding expansion coefficients 
are the expectation values Vij (/x) of the interaction potential V{r), defined by Eq. (^): 

N 

V{t) <pf\r) = Y: Vl[\f^ 0f (r) , £ = 0, 1 . (36) 

j=0 

Inserting these expansions involving the expectation values Vij (/i) of the interaction 
potential V{r) in Eq. (^, applying the Fourier-Bessel relation (^), and remembering 
the orthonormalization relation (^) satisfied by the momentum-space basis functions 
(f)i (p), the term Dij reduces to the — at least, for all power-law potentials algebraic — 
expression 

N N N 
r=0 s=0 t=0 
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B.6 Relations Among the Expansion Coefficients 

The expansion coefficients bij, Cij, dij defined by Eqs. (|28D, (|3TD , (|32|), respectively, are, 
of course, not independent; they satisfy the (only in the limit N -^ oo, exact) relations 

N 

/ . C^i Crj = Oij , (37) 

r=0 

N N 

H CU drjifi) = Yl ^h(/^) Crj = lififJ^) = bijifi) , (38) 

r=0 r=0 

AT 

^ci:,(/i)ci,,(/i)=/f(/i). (39) 

r=0 

Clearly, relation (^) expresses only the unitarity (in the limit A^ —>■ oo) of the matrix 
c = (cjj) which represents, according to its definition (^Tj), the transformation from the 
basis {0j- {p),i = 0, 1, 2, . . .} to the basis {(f)i {p),i = 0, 1, 2, . . .} of the space L2{R^). 

B.7 Accuracy Considerations 

The existence of relations like the ones in Eqs. (|37[) , (^8]), and (|39[) opens a possibility to 



investigate systematically the errors induced by the truncations of the expansion series 
in the terms Bij, Cij, and Dij. For instance, for d = 15 (i.e., 15 x 15 matrices) and for 
iV = 49 (i.e., a truncation to the first 50 basis states), the relative error is less than 3%. 
Furthermore, taking advantage of the fact that the momentum-space basis function 
0Q (p) may be cast into the simple form 

4°'(P) = 4,/^ ""' 



TT (p2 + yU^)^ 

the first entries of the terms Bij and Cij defined in Eqs. (^ and (0), respectively, can be 
evaluated, for appropriate choices of the interaction potential V{r), analytically. These 
explicit results may be compared with the corresponding approximate values obtained 
from the series expansions ( pOD and (|35|) , respectively. For the linear potential ([T0|) , the 
numerical value of the term i?oo is found to be 

oo oo 

^Eoo ^:^Jdkk'<j>^^\k)Jdk'k"Vo{k,k')<j>S'\k') = ^ = 21.3 . 



This has to be compared with the series expansion 

^ r=0 "^ r=0 

which yields, for, e.g., N = 49, that is, taking into account the first 50 basis functions, 

— E ^ro(l) v}o\l) = 21.333333333333334562 . 

Similarly, for a constant potential, i.e., V{r) = A, the numerical value of the term Cqo 
is found to be 

CO oo 

i^Coo = ^ fdke<pf{k) Uk'k''V,{k,k')<pf{k') = ^ = 21.3 . 
Au ir A u J J 3 

'^ '^0 
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This has to be compared with the series expansion 

g^ N N N 

■^ P r=0 s=0 r=0 

which yields, for, e.g., A^ = 49, that is, taking into account the first 50 basis functions, 

49 

Svr ^ clodroil) = 21.333307054 . 

r=0 

Beyond doubt, this accuracy of the series expansions should suffice for our purposes. 
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